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In the layer-by-layer controlled artificial structure of transition-metal oxides, it is known that 
the transition metal d charge valence, orbital occupation, and other degrees of freedom are often 
reconstructed, which can lead to unexpected material properties such as a highly metallic interface, 
magnetism, and superconductivity. In this study, by using first-principles density functional theory 
calculations for (LaNi03) m /(SrTi03) n superlattices, we report a different type of reconstruction 
caused by the heterostructuring, that is, redistribution. The Fermi level density of states of metallic 
nickelate layers is significantly reduced without charge transfer in the vicinity of interface to the 
insulating SrTi03. Further analysis indicates that this systematic way of electron response can 
be a general feature particularly in the nickelate systems. Our result sheds new light towards 
understanding the nickelates and other transition-metal oxide heterostructures. 
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I. INTRODUCTION 



Recent advances in the layer-by-layer growth tech- 
nique of transition-metal oxide (TMO) heterostructures 
have created considerable research interest [1]. In these 
artificially structured interfaces, many exotic material 
characteristics have been reported that are strikingly 
different from the bulk properties and from the char- 
acteristics of a typical semiconducting interface (non- 
TMO interface) [2]. For example, early studies of Ti- 
based superlattices such as LaTi0 3 /SrTi0 3 (LTO/STO) 
and SrTi0 3 /LaA10 3 (STO/LAO) have revealed that 
heterostructuring induces unique electronic processes in 
TMOs that can dramatically change the macroscopic ma- 
terial properties [J-Q. In these cases, it is widely believed 
that the valence charge reconstruction and the polar dis- 
continuity can drive the interfaces to be metallic even 
if the two mother materials are good insulators. The 
orbital degrees of freedom are also found to be recon- 
structed at the interface between manganite and cuprate 
[7|. Other examples of emergent material properties 
caused by the heterostructure geometry include super- 
conductivity Q and magnetism [3, LLQ| • These new find- 
ings raise the important question how the heterostructur- 
ing of TMO materials affects the density of states (DOS) 
close to the Fermi level and the related low-energy prop- 
erties. 

In this study, we report a new type of interface phe- 
nomenon: Fermi level DOS modulation caused by the 
state redistributions without valence change. Although 
this phenomenon exhibits similarities to the charge recon- 
struction and orbital reconstruction, it displays distinc- 
tive characteristics. We found that the Fermi level DOS 
strongly modulates as a function of the location of nickel 



'Electronic address: mj.han@kaist.ac.kr 



atoms relative to the interface. Our first-principles calcu- 
lations of LaNi0 3 (LNO)/STO clearly demonstrate that 
this novel electronic modulation around the Fermi level 
is intrinsic to the structure itself, i.e., the heterostruc- 
turing of nickelates, and not originating due to the other 
effects such as the charge transfer, oxidation of Ni and 
valence change. Further analysis strongly indicates that 
this kind of behavior can be a universal feature in nicke- 
late superlattices sandwiched by any wide-gap material. 
Our result sheds new light on the understanding of TMO 
interface phenomena. 

It should be emphasized that our system 
(LNO) m /(STO)„ is distinctive from a widely-studied 
system, LNO/LAO, due to the possibility that Ti can 
have active d electrons around the Fermi level. The 
charge transfer can in principle take place between Ti 
and Ni. As an example, it would be instructive to 
compare LNO/STO to LNO/LTO superlattice in which 
the Ti can clearly have an electron in its d orbitals, 
and there may be the electron transfer from Ti to Ni, 
leading to the configuration of Ti as d 0+s (or d x ~ s 
depending on the amount of charge transfer) and Ni 
as d 8 ~~ s (or d 7+s ). This kind of charge transfer is 
basically the mechanism that drives the exchange bias 
effect in the LNO/LaMn0 3 (LMO) [III- However, in 
our case LNO/STO, it is not at all clear if such an 
electron transfer would take place especially for the 
case of m=n=l. Note that (LNO)i/(STO)i can also 
be identified as (SrNi0 3 )i/(LTO)i, and this alternative 
specification of the system implies that Ti has d 1 
(LTO-like) configuration instead of d° (STO-like). As 
we will show in the below, this possibility is not realized 
and Ti remains as d° even in the m—n—1 case, and 
Ni as d 7 (we drop off the indication of ligand hole for 
simplicity). From this interesting new finding, important 
questions arise: If there is no charge transfer between 
the transition metal (TM) ions, what happens at the 
interface? In other words, under the condition that no 
charge transfer is allowed, what kind of response can 



be made by the metallic material, LNO, in the vicinity 
of interface to the insulating STO? The interface Ni 
would exhibit the same characteristics with the inner 
layer Ni? If not, what kind of possibility does Ni have in 
this superlattice geometry? This is a well defined open 
question that has never been addressed before. 

Our calculations show that the Ni-d electrons are actu- 
ally adjusted at the interface in a systematic way that its 
DOS is redistributed while keeping the same total num- 
ber of electrons. That is, the interface Ni, closer to the 
insulating STO, becomes more insulator-like in the sense 
that its Fermi level DOS gets reduced while the inner 
layer Ni more metal-like in the sense that the more DOS 
at the Fermi energy. And, we suggests, this DOS redistri- 
bution process can be understood from a simple picture 
based on the interactions between the heterostructuring- 
induced molecule-like states. 



Result 



Classification of Ni types and modulation of Fermi 
level DOS 



Let us start with a classification of nickel sites in 
LNO/STO considering their local environments; nickel 
can be located either in between two La 3+ 2_ layers or 
in between the La 3+ 2_ and Sr 4+ 2 ~ layers. Further, 
its neighboring TM can be either titanium or nickel (see 
the inset of Fig. Q}. If a nickel ion is located in between 
two La 3+ 2_ layers and its two neighboring TM sites are 
both Ni, it is bulk-like. However, at the interface, there 
are three different local structures. Let us denote these 
three interfacial types as II, 12, and 13 (see the inset of 
Fig. [1]). As the electronic structure of Ni can be affected 
by both these two factors (the neighboring TM ions and 
the ionic potentials caused by ^4-site cations), it is im- 
portant to understand these effects on the Ni electronic 
structures. We will discuss this point below in further 
detail. 

Recent depth-resolved standing-wave x-ray photoemis- 
sion spectroscopy (SWXPS) studies [12( show that the 
near-Fermi-level nickel states in (LNO)4/(STO)3 are sig- 
nificantly suppressed at the interface. A detailed analysis 
of the angular behavior of the SWXPS spectra indicates 
that the electronic states are suppressed at the two outer 
LNO layers adjacent to the STO layers (type 12 and 13 
nickel in our definition) , but not in the inner layers (bulk- 
like nickel). It was speculated in Ref. QJ] that the vari- 
ation in Fermi level DOS could be due to the different 
oxidation of Ni ions; however, no further experimental 
evidence or explanation was provided. 



DOS redistribution without charge transfer 

We note that a clear understanding of this DOS mod- 
ulation is of significant importance. If there is a charge 



transfer between TM ions, it is quite natural to ex- 
pect that the amount of DOS changes at the Fermi 
level. Interestingly, however, our calculation shows that 
such a charge transfer and oxidation process are not re- 
sponsible for this DOS modulation. We performed ex- 
tensive calculations for (LNO) m /(STO) n with (m,n) — 
(1,1), (2, 2), (3, 3), (4, 1), (4, 3), and (5, 2), and found that 
the number of e g electrons is 2.36-2.39 regardless of the 
type of nickel. From the point of view that all the nick- 
els remain in the same charge status, the observed DOS 
modulation is quite unexpected and hard to be explained 
by any of the 'reconstruction' processes that were sug- 
gested before. 

The fact that there is no valence change or further oxi- 
dation of Ni is most dramatically seen in m=n=l case for 
which the alternative characterization of the superlattice 
as (SNO)i/(LTO)i is also possible if the charge trans- 
fer might really happen. Since Ti has d 1 configuration in 
LTO while d° in STO, the result can be most clearly seen 
in the Ti states. The calculated Ti-d DOS in Fig. [2(b) 
shows that all the Ti-d states are empty. Therefore one 
can safely identify the system as (LNO) \j (STO) i, and 
there is no charge transfer between Ni and Ti even in the 
m=n=\ case. 

The same point is also confirmed by that the effect 
of A-site cation potential asymmetry (i.e., Sr 4+ versus 
La 3+ ) does not cause any significant change in the elec- 
tronic structure. This effect can be studied by examin- 
ing the (3,3) heterostructure, for example, that contains 
both type 12 and 13 nickel ions (compare the dotted-blue 
with dashed-green line in Fig. [3(b)). Since the two nickel 
ions only differ due to the ^4-site cations, the different 
electronic structure reflects the effect due to the ^4-site 
asymmetry. Figure [3(b) shows that the DOS change is 
very small. Therefore, the valence reconstruction is not 
relevant to this nickelate system presumably due to the 
strong covalency. In addition, for other (m,n) structures, 
the projected DOS for type 12 and 13 nickel is always 
similar (seen in Fig. H(b) for the case of (5,2) or upon 
comparison of Fermi level DOS value at z = and z = 1 
in Fig. H|). 

Now, considering that the nickel valence states are all 
same regardless of its position and the cation potential 
effect, it is quite surprising that the Fermi level DOS sig- 
nificantly varies. The calculated value of this quantity is 
summarized in Fig. Q] where the DOS values at the Fermi 
level (projected onto the different nickel sites) are plot- 
ted. It is clear that the interfacial nickel ions (see Fig. Q] 
where the z-position is normalized so that types II, 12, 
and 13 located either at z = or at z = I) always have 
notably smaller Fermi level DOS than the bulk-like nickel 
sites (0 < z < 1 in Fig. Q]) . The Fermi level DOS is sig- 
nificantly enhanced for the bulk-like nickel ions. The dif- 
ference becomes even more pronounced as the number of 
Ni layers increases from m = 3 to m = 5. For m = 5, the 
Fermi level DOS for the inner most layer Ni ions is larger 
by a factor ~2 when compared with that of the interfa- 
cial one. Since there is no charge transfer and no valence 



change, but a strong DOS modulation at the Fermi level, 
DOS should be redistributed in such a way that the total 
number of electrons is kept same. The further analysis 
shows that DOS is actually redistributed in a systematic 
way so that the Fermi level DOS gradually decreases as 
the nickelate layer gets closer to the interface. In order to 
make the charge valence unchanged, this weight is trans- 
ferred to the lower energy part as schematically shown in 
Fig. SJa). It is interesting to note that the metallic LNO 
layers become more insulator-like when they come closer 
to the insulating STO, while keeping the same valence 
charge; i.e., they try to have smaller number of states 
at the Fermi level. On the other hand, the inner LNO 
layers, far away from STO, become more metallic in the 
sense that they have more DOS at the Fermi level. 

This process is schematically illustrated in Fig. EJa), 
and Fig. 0Jb) shows that the reduction of DOS is quite 
significant. As seen in Fig. 2Kb), the d^-r 2 states of 
the inner most (B; bulk- like) nickel is much larger than 
the two interface nickels (12 and 13). The arrow indi- 
cates the DOS reduction at the Fermi level. It is found 
that the reduced states are transfered to the lower energy 
parts in case of interface nickel sites (for the wider range 
plot and the integrated DOS are given in Supplementary 
Material). 



Mechanism of redistribution 

It now became clear that the Fermi level DOS behaves 
in a systematic way of state redistributions while keeping 
the same total number of charge. The remaining question 
is how to understand (at least qualitatively) the mecha- 
nism behind the DOS modulation or redistribution pro- 
cess. We suggest a relatively simple picture can explain 
this novel interface phenomenon based on the interac- 
tions between the molecule-like quantum states formed 
by the heterostructure geometry. Let us go back to the 
m= n=l case which serves as a starting point of our ar- 
gument. The Ni-e 9 DOS for to = n = 1 superlattice 
is presented in Fig. Ufa). The important point here is 
that the bandwidth of d 3z 2_ r 2 state is markedly narrower 
than that of the d x i_ y 2 state. The smaller dispersion of 
d 3z 2_ r 2 is due to the limited hybridization along the z 
direction by the presence of STO layers. From this char- 



acteristics of d 32 2_ r 2 state, we can consider this state as 
remnants of the molecular orbital character arising from a 
chain of to d 32 2_ r 2 orbitals along the z direction. These 
molecular-orbital-like features are further confirmed by 
the DOS shape of the (3,3) case; Figure [2Kb) clearly 
shows the features related to the bonding, antibonding, 
and nonbonding states of the to = 3 chain. It is noted 
that the DOS of the bulk-type nickel (solid-red line) is 
strongest in the bonding and antibonding features at — 
and -2.2 eV. 

This point is further supported by our tight-binding 
analysis, which captures the essence of the features found 
in the first-principles results. Assuming that the sys- 



tem can approximately be treated as a linear chain of m- 
molecular orbitals, we constructed a tight-binding model 
Hamiltonian corresponding to N atoms, H = J2i=i t c l c j> 
where c\ creates an electron at the ith site in the z direc- 
tion. The weight of the central (bulk-like) atom |(i c |i?}| 2 
is plotted in Fig. [3ja) as a function of energy E (in the 
unit of t). As the chain length, to, increases from to = 1 
to 3 and 5, the states interact with each other and spread 
out over a wide energy range. For m = 3, the two-peak 
feature of bonding and anti-bonding is most clearly seen, 
which is consistent with the results of our first-principles 
calculation in Fig. [3Kb) (for the results of other to val- 
ues, see Supplementary Materials). Another important 
point in this analysis is that the bonding-peak position 
gradually shifts toward the lower energy region as to in- 
creases, as marked by the asterisks in Fig. EK a )- Figure 
[3Kc) clearly shows that this characteristic behavior of the 
lower energy part of the d 3z 2_ r 2 states is also found in 
the first-principles calculations. As the number of lay- 
ers increases, the peak position (marked by asterisks in 
Fig. [3Kc)) shifts towards the lower energy region. It is 
the electronic origin that leads to a strong modulation in 
the number of states at the chemical potential shown in 
Fig.[TJ Therefore, we can conclude that the characteristic 
modulation of Fermi level DOS is attributed to the for- 
mation of molecule-like quantum states induced by the 
blocking of hybridization along z direction, and to the 
redistribution of DOS around the Fermi level caused by 
the interactions between these states. 



Discussion 

There are crucial differences between the DOS modu- 
lation and the previously reported interface phenomena. 
In general, electronic Q and orbital [3] reconstructions 
depend on the relative on-site energies of interfacial ions 
leading to a redistribution of charge between different 
sites and orbitals. In contrast, Fermi-level DOS mod- 
ulations involve a strong redistribution of states within 
the same sites and the same orbitals due to the presence 
of the interface. The redistribution leads to a strong 
change in the contribution to the Fermi level DOS as 
a function of the position of the TM relative to the in- 
terface. Since the low-energy properties are mainly de- 
termined by the states close to the chemical potential, 
strong position-dependent electronic properties can be 
expected and macroscopic properties (e.g., resistivity) 
should be sensitive to these strong modulations. The 
direct observation of this kind of electronic response is 
difficult as the charge valence of the interface remains 
unchanged, and only the near-Fermi-level DOS should be 
measured. However, we note that SWXPS measurements 
can be used for the verification of our results. Therefore, 
it is important to note that the results of a recent experi- 
ment by Kaiser and co-workers for (m=4, n=3) compares 
well with our conclusions [12j , 

Our analysis can be applied to other superlattice struc- 



tures with LNO sandwiched by any wide-gap insula- 
tor in which the hybridization is also strongly blocked 
only along the z direction. Therefore, the same type 
of DOS modulation should be found in the related het- 
erostructures. To verify this point, we calculated the 
Fermi level DOS for (LNO) m /(LAO)„ with (m,n) = 
(3, 3), (4, 1), (5, 2), and we found that the same modula- 
tions also occur in this different type of nickelate systems 
(for the plot comparable with Fig. [1] see Supplementary 
Material). These observations further confirm our anal- 
ysis and conclusion. It is noteworthy that our tight- 
binding analysis based on the molecular-orbital chain 
model is better applied to the case of LNO/LAO. We 
can conclude that the DOS modulation can be a univer- 
sal feature in nickelate superlattices and possibly also in 
other interfaces in between metals and insulators. 

In addition, our results concerning the Fermi level DOS 
reduction at the interfaces can provide an interesting new 
picture of nickelate superlattices. It is noteworthy that 
the DOS reduction at the interface is observed for all the 
compositions of (m >3, n), while for m < 2, only the in- 
terface nickel ions exist and there is no bulk-like one. In- 
terestingly, several experiments on the nickelate superlat- 
tice independently reports that the MIT occur at to w 3. 
The same critical thickness was commonly observed in 
the nickelate systems; LNO/STO (Tl, LNO/LAO Q, 
LNO/SrMn0 3 0, and LNO/LMOlll. It is therefore 
temping to relate the metallic and insulating phase to the 
bulk-like (B) and interface-like (II, 12, 13) nickel layers in 
our classification, respectively, as the bulk-like nickel has 
more states at the Fermi energy and it start to appear 
at to = 3. Further, a recent experiment by Boris and 
co-workers reported that this MIT is accompanied with 
the paramagnetic to magnetic transition [14(. We note 
that, in the Kondo-type screening, the screening strength 
is governed by the Fermi level DOS as is clear from 
the Kondo temperature scale, Tr- ~ € ~ 1 /poJk ^ wnere p Q 
refers to the Fermi level DOS. Therefore, it may be in- 
teresting to regard the nickelate superlattices based on 
our result as a kind of 'Kondo lattice' system where the 
enhanced Fermi level DOS is responsible for the metallic 
conduction and simultaneously for the screening of local 
magnetic moments (schematic illustration of this point 
can be found in Supplementary Material). 

To summarize, we report a new type of interface phe- 
nomenon: Fermi level DOS modulation across the LNO 
layers in LNO/STO superlattice systems. We show that 
this modulation is an intrinsic effect caused by the het- 
erostructuring itself without valence change or oxidation. 
It is related to but clearly distinctive from other intrigu- 
ing interface phenomena such as orbital and charge re- 
construction. Our analysis demonstrates that this novel 
electronic response to the hetero-interface structure orig- 
inates from the novel process of quantum state formation 
and interactions between them. It can be responsible for 
MIT in the related systems, providing a new theoretical 
aspect to the complex-TMO research. 
Methods 



rection, and a linear combination of localized pseudo- 
atomic orbitals (LCPAO) as a basis set. We adopted 
the local density approximation (LDA) for the exchange- 
correlation energy functional as parametrized by Perdew 
and Zunger [17|, and we used an energy cutoff of 400 
Ry and k-grid of 12 x 12 x 6 per unit volume. It is 
noted that bulk LNO is a paramagnetic metal down to 
low temperatures and no studies yet have reported on 
the magnetic order in LNO/STO. Therefore, a conven- 
tional LDA calculation is appropriate for the LNO/STO 
system since a LDA+U type of calculation, for example, 
has to assume a long-range ordered magnetic phase. We 
used the Mulliken method for the charge analysis where 
the Kohn-Sham states are projected onto our LCPAO 
basis orbitals. It is noted that the number of Ni-e ff elec- 
trons obtained in this study is in the range comparable 
to the case of LNO/LXO (X: B Al, Ga, In) superlattice 
systems as reported in Ref . [l8| [l8| . The geometry relax- 
ation has been performed with the force criterion of 10 -3 
Hartree/Bohr. During the relaxation process, we assume 
that the in-plane lattice constant does not change due 
to its pinning with the substrate (we assumed the STO 
lattice parameter for the substrate) [19j . 
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For band structure calculations, we employed the 
norm-conserving pseudopotential with a partial core cor- 
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Figure Captions 

1. The calculated Fermi level DOS for various (m,n) 
combinations of superlattices. The z-axis positions 
of Ni ions are renormalized so that the first Ni0 2 
layer is located at z — and the last is located at 
2 = 1. (Inset) Four different types of Ni sites exist- 
ing in the (LNO) m /(STO) n superlattice. Dark gray 
and light blue circles represent Ni and Ti, respec- 
tively. Solid and dotted vertical bars in between 
transition metals represent SrO and LaO layers, re- 
spectively. Nickel of type II is located in between 
the Ti0 2 -SrO and LaO-Ti0 2 layers, type 12 in be- 
tween Ti0 2 -SrO and LaO-Ni0 2 , and type 13 in be- 
tween Ni0 2 -LaO and LaO-Ti0 2 . Nickel of type B 
is in the bulk-like arrangement located in between 
Ni0 2 -LaO and LaO-Ni0 2 . 

2. (a) Projected Ni-e g DOS in the (1, 1) superlattice. 
Red (solid) and blue (dashed) lines represent the 
d^ z 2_ r i and d x 2_ y 2 states, respectively, (b) Pro- 
jected Ti-£2g DOS in the (1, 1) superlattice. Blue- 
dashed and red-solid lines represent the d xy and 
d yz ,zx states, respectively. Vertical dotted lines in- 
dicate the position of Fermi energy. 

3. (a) Tight-binding calculation results for a linear 
one-dimensional molecular chain for the weight of 
the central atom |(i c |.E}| 2 in the energy eigenfunc- 
tions \E). The length of the chain corresponds 
to the thickness of the LNO layer, i.e., m in 
(LNO) m /(STO)„. The central position is given by 
i c = rn/2 and (m + l)/2 for even and odd m, re- 
spectively. Energy unit is in the hopping param- 
eter t. (b) Projected d 3z 2_ r 2 DOS for (3,3) case. 
As expected, the bonding, non-bonding, and anti- 
bonding characteristics are most clearly seen in the 
case of (3,3). Type 12, 13, and B states are repre- 
sented by blue-dotted, green-dashed, and red-solid 
lines, respectively. The Fermi level is set to be zero. 
(c) The projected d 3z 2_ r 2 DOS for bulk- type Ni in 
the (3,3), (4,1), and (5,2) structures. For compari- 
son, type 12 in (2,2) is also shown (red). The Fermi 
level is set to be zero. The asterisk marks indicate 
the peak positions. 

4. (a) Schematic illustration of DOS reduction at the 
Fermi level as a function of layer position relative 
to the interface. The reduced DOS is transfered to 
the lower energy part, (b) Projected d 3z 2_ r 2 DOS 
of two interface layer nickels (blue and green) and 
the inner layer nickel (red) in the (5,2) superlattice. 
The arrow indicates the amount of DOS reduction. 



O 
Q 

"oi 

> 



0.8 



0.7 



0.6 



0.5 



0.4 




0.2 0.4 0.6 0.8 

Normalized z-position of Ni 

FIG. 1: 



Fig.l 





U.fa 


1 1 r 


i 2 _ j> i 






0.4 


. (a) Ni 




3z 2 -r 2 


eg 


> 

CD 










A 


CO 
(1) 


0.3 








■H 
V 














m 
























</> 














O.J' 








U) 










O 






I ». ••' I '• 




Q 


0.1 





















> 








cu 








u> 




1 






a> 






















: X 




CO 




y 


CO 


> 

CD 


0.8 


(b) Ti 






V z ' 


co 

o 


CO 
CD 


0.6 








•II 
















ffl 








V"/ ■ 1 




-t— • 








V : 








0.4 












U) 






. 








O 














> 


a 


0.2 












-9? 

a> 
to 









, •" , J 




; i j 


CO 




4 -2 


2 4 § 








En 


ergy (e\ 







□ 



FIG. 2: 



0.5 



(a) 


1 m=1 


* 
* 





-2t 



2t 




Fig.2 



Energy (eV) 
FIG. 3: 



Fig.3 



(a) 



Fig.4 



reduction 



Inner layer 




Fermi Level 



(b) 




1 2 

Energy (eV) 



FIG. 4: 



